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TECHNICAL NOTE 4026 

EFFECTS ■ OF THERMAL RELAXATION AND SPECIFIC-HEAT CHANGES 
ON MEASUREMENTS WITH A PNEUMATIC-PROBE PYROMETER 
By P. W. Kuhns 


SUMMARY 

Equations are derived which evaluate the effects of thermal vibra- 
tional relaxation and specific-heat changes on pneumatic -probe pyrometer 
measurements in the region 1000° to 3000° K and 0.3 < Mach number <2.0. 
Examples are given for typical probes in two combustion mixtures. 


INTRODUCTION 

During recent years much interest has been shown in the problem of 
the effects of thermal relaxation or disequilibrium of vibrational modes 
of excitation on various aerodynamic processes. The effect of the vibra- 
tional relaxation time on the thickness of shock waves has been investi- 
gated in reference 1, and the shock-wave thickness has been used to 
measure relaxation times in references 2 to 4. The effect of vibrational 
relaxation time on flow through a rocket nozzle is studied in reference 
5. The manner in which vibrational relaxation affects the measurement 
of total pressure has been used to find the vibrational relaxation times 
of some gases in references 6 to S. A discussion of the effects of re- 
laxation on same nozzle processes Is found in reference 10. 

A review of vibrational relaxation measurements is presented in 
reference 11 with a fairly complete bibliography. This report also con- 
tains a table of roam temperature values of the average number of colli- 
sions necessary to achieve vibrational relaxation as found by many in- 
vestigators. Since that time (l95l) , work has been done both theoretically 
(refs. 12 to 14) and experimentally (refs. 2 , 3, 4, 8, 9, 15, 16, and 17) 
in the determination of vibrational relaxation times. A more recent re- 
view of relaxation phenomena and their measurement may be found in refer- 
ence 18. 

This report considers the problem of the effects of vibrational re- 
laxation on pneumatic -probe temperature measurements. The pneumatic- 
probe pyrometer does not directly measure the gas temperature but measures 
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other aerodynamic quantities and relates them to the gas temperature by 
means of theoretical equations. The basic formulation of the method may 
be found in references 19 to 23. Some experimental results using the 
method are reported in reference 21. Reference 23 presents data compar- 
ing the results of the method with simultaneous data from thermocouple 
and sodium D line reversal measurements. 

As the gas flows from the free stream into the pyrometer., it suffers 
a partial stagnation because of the presence of the probe body and an 
acceleration because of the pressure drop in the probe. The result is a 
rapid change in the available thermal energy. The relaxation (response) 
times of some vibrational modes of thermal excitation are such that the 
modes cannot follow these changes. As a result , the critical velocity 
in the probe and the temperature as calculated are affected. 

Equations are derived herein which can be used to correct the meas- 
ured parameters in such a way as to obtain the proper temperature when 
the vibrational modes of a molecular species are only partially relaxed. 
By using these equations, the effects of thermal relaxation and specific- 
heat changes on the measured total pressure, the calculated total temper- 
ature, the critical velocity, and the mass -flow rate are considered. 

The equations are evolved in a general form so that they may be appli- 
cable to similar flow problems. 


THEORY 

Outline of Subscripts and Superscripts Used 

As an aid in reading the text, the following guide to the subscripts 
and superscripts is presented. A complete list of symbols is presented 
in appendix A. 

Subscripts . - The magnitude of many of the quantities used is de- 
pendent upon the state of relaxation or equilibrium of the modes of 
thermal excitation. For convenience it is possible to group the modes 
of thermal excitation of all the different molecules in a gas according 
to their probable state of thermal equilibrium or disequilibrium at the 
first nozzle of the pneumatic probe (fig. l) . 

The first subscript denotes the relaxation group into which the 
quantity is to be placed. 

none all modes 


P 


relaxed modes 
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There are extremes in the quantities grouped as relaxed modes. 

These extremes occur when the partially relaxed vibrational mode approaches 
either full relaxation or an unrelaxed state. These extremes are denoted 
by: 

R relaxed modes when partially relaxed mode is assumed relaxed 

U relaxed modes when partially relaxed mode is as sinned unrelaxed or 

"frozen” 

Z relaxed modes when partially relaxed mode is assumed unrelaxed and 
specific heats of all modes are assumed to be independent of 
temperature 

The other modes are designated as follows: 

B partially relaxed vibrational modes 

X vibrational modes which are considered unrelaxed or frozen 
Thus , 


E = Ep + Eg + Ej^ 

A second subscript denoting the station where the quantity is to be 
taken is shown as 0, 1, 2, 3, or 4. 

Superscripts . - Usually the state of relaxation also denotes the 
temperature at which a quantity is considered. When there are exceptions, 
the temperature state at which the quantity is to be taken is given by a 
superscript . 

Thus, Eb is the energy of the B mode taken at the temperature of 
the B mode, while Eg is the energy of the B mode taken at the tem- 
perature of the P modes. 

There is also a difference between the measured quantities and the 
quantities involved in the theoretical equations. The measured quanti- 
ties will be designated by the superscript *. 


Theory of Operation 

In the pneumatic -probe method of gas temperature measurement, the 
total temperature is usually found by using the following procedure, 
with reference to figure 1: 
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(1) A pressure drop through the nozzles is established by an exhaust 
system large enough that the flow through both nozzles becomes critical. 
The mass flow reaches a maximum value (ref. 23). 

(2) After passing through the first nozzle, the hot gas is cooled 
before reaching a second nozzle. 

(3) The total pressure P4 and the total temperature T* at the 
second nozzle are measured. 

(4) The flow in the probe is halted and the total pressure Pg is 
measured at the first nozzle using the probe as a pitot tube. 

(5) The true total temperature Tq of the gas is obtained as a 

function of the above measurements from aerodynamic relations equating 
the mass flows through the two nozzles. 


Velocity Profiles, Relaxation Times, and Specific-Heat Changes 

In making the measurements outlined in the preceding section the 
gas flow is changed somewhat. Such changes in gas flow can, under some 
circumstances, affect the measurements. 

As the gas passes from the free stream into the probe, it is first 
slowed down by the presence of the probe and then speeded up by the pres- 
sure drop at the nozzle. This is shown in figure 2, which presents 
measured velocity profiles below a free-stream Mach number Mp = 1 for 
the center streamline of a simulated probe and shows the probe shape 
drawn to scale. 

The inequality in the decelerating and accelerating distances (and, 
thus, times) should be noted in figure 2. For the typical probe size 
(D ■ l/2 in.; d * l/8 in.) in a gas at 2500° K with a free-stream Mach 
number near 1, the decelerating time is about 10 microseconds; the ac- 
celerating time Is about 1 microsecond. 

If the components of a gas are to remain in equilibrium during this 
process, energy must be exchanged between molecules and between modes of 
thermal excitation of the molecules. If this exchange takes longer than 
the accelerating time, the gas will not be In complete equilibrium, and 
the flow through the nozzle will be affected. 

When the pressure P* is being measured, the flow through the 

nozzle is stopped. The velocity profile now is like that shown In fig- 
ure 2 for the velocity ratio Ug/u^ =0. If the time a gas molecule 
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takes to stop in front of the probe is on the same order of magnitude 
as the time to reach equilibrium, the measured pressure Pg will be 
affected. 

At the temperatures, pressures, and velocities to which the probes 
are usually subjected, the rotational and translational modes of thermal 
excitation establish equilibrium in less than 1/10 microsecond so that 
they can follow the energy changes. However, some vibrational modes 
cannot follow the changes, having relaxation times t of 0.1 to 
10,000 microseconds. 

The effect of this relaxation time upon temperature is shown in 
figure 3, in which a gas mixture is assumed that, at the first nozzle, 
has in its constituents one of three vibrational modes: one that is re- 

laxed (in equilibrium), a vibrational mode that is unrelaxed (unaffected) 
and one that is partially relaxed. The temperature and velocity at the 
second nozzle are usually so low that any effects due to relaxation may 
be neglected. 

It should be noted in figure 3 that not only doeB the relaxation of 
a vibrational mode affect the temperature of that mode, but, because 
energy must be conserved, it affects the temperatures of all relaxed 
modes . 

The time to reach equilibrium is measured by the relaxation time, 
which is defined by the equation (refs. 11 and 18): 

®B “ “ ^ CD 

for constant vibrational specific heats. 

To write equation (l) it was assumed that the specific heats were 
constant during the energy change. However, by using a very rough calcu- 
lation it is easily seen that even at subsonic speeds the temperature 
of the gas may vary as much as 20 percent when the gas passes from the 
free stream into the probe. In the temperature range from 500° to 3000° 
K the vibrational specific heat of most gases will change with tempera- 
ture. This is shown in table I which lists values obtained from refer- 
ence 24. For this reason equation (l) should be modified to: 

3E % - - £ 0% - 4 ) (2) 

where is the energy of the partially relaxed vibrational mode taken 

at the temperature of. the relaxed modes. 
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Derivation of Equations for Tq 

The true total temperature Tq is defined as the temperature the 

gas would have if it were brought adiabatically to a stop with all modes 
in equilibrium: 

2 


Et + 


U n 


t 0 s e 0 = 


4 /" 


( 3 ) 


This is the temperature which the probe is presumed to measure. 

As an aid in calculation, a second temperature will be defined: 

2 


u. 


T Z,1 = S 1 + 2? 


--eJ 

>,1 ' 


1 + 


'P,l 


- 1 


) 


U) 


where 


y t 2 x,2 1 

and T z ^x is not the total temperature but in many cases is close to it. 

In the derivation of the relation between Tq and the measured 
quantities P*, P4, and use is made of five fundamental equations 

which must be satisfied. — 

(l) The equation of conservation of energy: 


u 


u=r 


E B + % + %,1 + ~2 m ®B,1 + %>,1 + ®X,1 + — = C P,1 T Z,1 + ®B,1 + ^,1 

( 5 ) 

(2) The equation for critical flow: 


dm 5 




dpg ° = &P4 

(3} The equation for entropy change: 


** -/"* (4 ' i 


( 6 ) 


(?) 
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(4) The relaxation equation (eq. (2)) 

(5) The equation for conservation of mass: 

m 4 = m 2 = A 2 p 2 u 2 =* A 4 p 4 u 4 (8) 


In order to find a usable solution to these equations, four assump- 
tions are made: 


(l) The specific heats vary linearly in the form 


c i = c i,l + 2c i,l*i 



( 9 ) 


where is a constant and i signifies any particular mode. 

( 2 ) The static pressure at station 2 is related to that at station 1 
by the function 


p i Iv 

(3) The relaxation time t is constant over the process in front 
of the nozzle. 

(4) At station 1 all modes are in thermal equilibrium: 

9 1 “ ®P,1 = ®B,1 = ®X,1 C 11 ) 

By using the previous five basic equations and the four assumptions, 
an equation for Tq is found: 

/a 2 p |\ 2 ^ ) 

T ° = fof) ^ ~^uT (1 + ' a)(1 + ' b!(1 + 6o) (12) 

where 

.(r) - r 


\n-l 


(10) 


(13) 
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and 

1 + e a correction to mass-flow rate due to relaxation and specific- 
heat change with gas passing through a nozzle 

1 + c-jj correction to total pressure due to relaxation and specific- 
heat change when the flow through the nozzle is stopped and 
Pg is measured 

1 + *c T o/ T Z,l 


Appendix B gives a detailed derivation of the equations for 1 + « a and 

1 + e-jj. In part of appendix B the additional assumption is made that 

.the equation for critical flow (eq. (6)) for the first nozzle is satis- 
fied hy a critical velocity given as 



^ T Z,1^ 



(14) 


Working Approximation of 1 + 1 + and 1 + e c 


Even when making an assumption as to the velocity (eq. (l4 ) ) , the 
computation of 1 + e a and 1 + hy theBe methods is still too complex 
to he of any use to the engineer. A working approximation was found 
that appears adequately accurate. 


When computing e a and e-^ from the relaxation equations, the 
variables Z a , 2^, T a; and appear. These variables are derived in 

appendix B and are given hy 



The velocity profile u/u-j_ used to calculate 2 & and is that 

with flow through the nozzle. The velocity profile used to calculate 
and r b is that with the flow in the nozzle stopped. 
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With assumed values for e x and Mg L the specific heats and 

values of t are computed from tabulated data (ref. 24) or from tables 
I and II. The value of Ug/u^ is computed from, equation (13) . Then 
® a is calculated from the form ula ; 


* * e 

a a 


5 s "- i] *$)* ^|rr 


(17) 


where ® a ^u and ® a ^ R are the two extremes of e &} and is com- 

puted from the velocity profiles (eq. (l4)). The functions * a R and 
e ajU can calculated from the following equations: 


where 


= F)1 AjjA ^ 


a > u " 


- 1 


1 + -£i 


Tp T ~ 1 


„ .. 

%2 - 29 i ? — 

T P,i - 1 - 


(18) 


(19) 


where 


R % “ ^l+P ^ 

°U,2 - ®1 0 p )1+E% + % 

2 

\7"l5 -I -1 


a,R [ Q r 


, 0 Z,2/ V 6 1 / 


(20) 


( 21 ) 


( 22 ) 


E % “%1 + (=P,1*E + %1*B> - l) ^ - |V CT1 ; (23) 
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„ = e. 


'P,i l e 


Z ^2 . ! 


+ 6, 


R, 2 3_ cp 2 _ RPj^ 

The curve for is given by 

®b = ®b,U + (1 - ^ - r b)(®b,R “ 6 b,u) 


(24) 


(25) 


and 


where 


and 


where 


2 e 


Pz .1 


- [%±\ R /_fl_\ 2PU . n 
1 %^/ ( ■%, 2 ) 


(26) 


T Z,1 = + M Z,l) 




“ C P,1*P l "If - * 


(27) 

(28) 


T, 


■Z,l 


'P,i\ a 


i± - 1 


^■ Sl =p, 1+ EPff 


+ Si 


(29) 


2 c 




v23? 


'b ; R \T e ^ 2 


V\ R {JkX 1 * _ , 


% 2 , 


(30) 


R{3 E “ C B,1 + + C B,1^B^ (^-ff " X j 


(31) 
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(32) 


By using the assumption expressed in equation (9), a relation for 
1 + £ c may be written in the following manner: 


1+ 6 c 


Of 


^ f 1 - jd 

C 1 + - i) 


C 1 


T, 


ZA 


(33) 


By using the above equations comparisons were made with the more com- 
plex formulas outlined in appendix B. The above equations are shown in 
appendix B to be accurate in the region 0.3 < Mp < 1.4. Modifications 

of the above formulas for the region 1.4 < M-j_ < 2.0 are given in appendix 
B (eq. (B58) ) . Below Mp = 0.3, Ug/up cannot be assumed constant. 


Other Corrections 

Specific-heat ratio . - If in the computation of hCrp^) some 
other value of r was used instead of Vp a correction can be made of 


1 


+ 


P! 


hCrpi) 

h(r) 


(34) 


Discharge coefficient . - In all the previous derivations the effec- 
tive areas A were used. These areas are equal to the discharge coef- 
ficient multiplied by Jtd^/4. The discharge coefficient Is affected by 
thermal relaxation. By using reference 23, an estimate of the error in- 
duced by neglecting this change was calculated. The maximum error varied 
frcan 0.04 percent at = 0 to -0.2 percent at Mp = 2.0. Thus, this 

effect may be neglected over the range of velocities and temperatures 
considered in this report. 

Dissociation . - In all the previous derivations and equations it 
has been assumed that the specific gas constant B was the same at both 
nozzles. Above 2000° K, dissociation of the gas products of combustion 
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make this assumption false. If differences in R are allowed., equa- 
tion (12) should he rewritten as 


T 


0 





" h ( r P,l) 

L"wtrJ 


(l + S a ) (1 + 8- b )(l + s c )(l + ® e ) 


(35) 


where . 


1 + e £ 



(36) 


Values of « e for three combustion processes at normal atmospheric 

pressure are shown in figure 4. The values in this figure were calcu- 
lated from reference 25 assuming that Q 4 < 1500° K. Below 2000° K, 

dissociation effects are negligible at atmospheric pressure. The quanti- 
ty 1 + ® e may be calculated for other pressures and gas constituents by 

using the dissociation constants found in references 25 to 27 or may be 
obtained graphically in the range of pressures of l/4 to 100 atmospheres 
by using reference 26. The composition of air at various pressures and 
temperatures up to 24,000° K can be found in reference 28. 

The effects of dissociation upon the specific heats and the entropy 
include not only the energies of each constituent but also include the 
energy of dissociation. The effect of such inclusions upon the thermo- 
dynamic properties of a diatomic gas are given in reference 29. Such 
effects must be included in the determination of e a and etj. 


The next question that arises at present cannot be answered with 
any degree of certainty. This concerns the dissociation relaxation 
times of the molecules. If dissociation is thought of as a large an- 
harmonic vibration, it would be expected that the dissociation relaxa- 
tion times would be comparable to the vibrational relaxation times at 
temperatures where there is appreciable (e.g., 10 percent) dissociation. 
The deviation of the vibrational relaxation time of oxygen from the 
theoretical curve at high temperatures may be given as evidence of this 
(fig. 5). However, this is contradicted by theoretical values obtained 
by using references 30 and 31. 

For this reason the region above 2500° K is usually avoided in this 
report. All experimental data concerning dissociation relaxation should 
be investigated before making any calculations in a temperature region 
of appreciable dissociation. 

v. 

Other methods of measuring Pj>. “ ® ie previous equations are de- 
rived with the assumption that ' is measured by stopping the flow In 
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the pneumatic prole . If the measurement is made by another probe in the 
moving gas stream, the equations for 1 + e-^ can still be used with the 

values of 2-^ and for the probe used in the equations. However, 

if P* is measured in the plenum chamber before a supersonic nozzle 
(P* = P Q ) , and, since in the derivation P^ was taken behind a super- 


sonic shock, an additional correction to 
for the entropy rise through, the shock: 



must be made to account 



APPLICATION 

Determination of Relaxation Times 

In order to use the previous formulas in determining the total tem- 
perature, it is necessary to find the relaxation times of the various 
components. The relaxation time depends upon the temperature, pressure, 
and the species of molecules colliding. The relaxation time of various 
molecules when colliding with another molecular species with a partial 
pressure of 1 atmosphere is shown in figure 5 as a function of the tem- 
perature. The data shown were taken from references 3, 4, 6, 7, 11, and 
18. Some of these data were originally published in the form shown, 
other data were converted from the number of collisions required to 
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obtain equilibrium. The conversion was made by using the relation 
t = (number of collisions }/(collision rate) and by computing the collision 
rate for collisions between hard spheres (ref. 33). Data on oxygen- 
nitrogen relaxation (not shown) may be found in references 4 and 11. 


In extrapolating the data for collisions between molecules of the 
same species , use was made of the relation first presented by Landau 
and Teller (ref. 34) and expanded by Bethe and Teller (ref. l) : 

, aSp-l/3 

t i„j “ e (40) 

(shown by solid lines in fig. 5) where 4 is the relaxation time of 
molecule i due to collisions with molecule j and is the partial 

pressure of molecule j . 


It has been found experimentally that for collisions between two 
different species the number of necessary collisions is usually almost 
independent of temperature. For these cases extrapolation can be made 
using the perfect gas law (ref. 33): 


T 






(41) 


(shown by dashed lines in fig. 5). 

In a gas mixture it may also happen that the vibrational modes of 
a molecular species are relaxed by more than one type of collision. For 
example , as seen in figure 5, the relaxation of nitrogen by both nitro- 
gen and water must be considered. In such cases the relaxation time is 
given by 



(42) 


as the collision processes are parallel. 

Figure 5 shows some disagreement between the theory (eq. (41)) and 
the experimental results for the relaxation of carbon dioxide by water. 
Much of this can probably be attributed to experimental error as relaxa- 
tion times were obtained from differences in relaxation. The deviation 
of the data for oxygen relaxed by oxygen from theory (eq. (40)) may be 
attributed to anharmonic vibration or dissociation. 

In all the above derivations and formulas it has been assumed that 
the relaxation time is constant. In fact, however, the relaxation time 
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Is dependent on the temperature 0p and the static pressure p. Use of 

such a variable relaxation time,, however, would make the computations 
more tedious than they already are. Below Mach 1 and at temperatures 
above 1500° K the free-stream temperature and pressure will give a fairly 
accurate relaxation time. This is true because the pressure and temper- 
ature will not vary appreciably below a free-stream Mach number of 1, 
and in the temperature range of interest, the value of t does not vary 
greatly. 

At supersonic speeds this simplification will not be valid. Here, 
all the relaxation effects take place after the shock wave; when the 
velocity is once again below Mach 1, pressure and temperature will again 
vary inappreciably. Thus, for speeds above Mach 1, the temperature and 
static pressure which should be used to determine the relaxation time are 
those of the stream after the bow wave. Denoted by the subscript T, 
the pressure and temperature to be used in determining t (from ref. 32) 
are given by 


9 ^. * 0-^, M-j_ < 1 


fl T" ®1 


[sr^ - (r^ - 1)] [(r x - i)Mf + 
(ri + 


, M. > 1 


(43) 


P T “ P l' **1 < X 
2y,M? - (r, - 1) 

Pt - Pi ■' Ml > X 


(44) 


Choice of Specific Heats 

After having determined the relaxation times of the various vibra- 
tional modes, it is possible to determine if such modes are to be con- 
sidered unrelaxed, partially relaxed, or totally relaxed in the previous 
equations for Tq. 

Whether the vibrational states of a molecule are relaxed or un- 
relaxed depends upon whether the relaxation time T is short or long 
compared with the time taken for the molecule to change energy (t f , fig. 

3) . This fact can also be written by using a relaxation parameter 
D/u^tt + which is approximately the ratio of these two times t^/v. Be- 
low D/uj_f + = 0.01 the vibrational energy of the molecule is "frozen" 



16 


MCA OT 4026 


at the initial -value, and e & and 6-^ are equal to 6 a ^u and u» 
respectively. Above D/uqt+ = 100 the vibrational energy of the mole- 
cule is in equilibrium with the other modes, and s a and e-g have the 
values s a ^R and Sq, g. Thus, when considering a gas mixture, the vi- 
brational specific heat of those constituents which have relaxation times 
such that D/u£t + < 0.01 will be given the designation c^, and these 
vibrational specific heats will not be used in the computation of c-g, 
cp, T 2 , 1 + e a , or 1 + 8-g. However, they will be included in c when 
computing Tq and 1 + e c . 

When a relaxation time is such that D/uqT + > 100, these vibrational 

modes will be considered fully relaxed, and the vibrational specific 
heats of their modes are included in the specific heats cp and c in 
all calculations. 

Vibrational specific heats with relaxation times such that 0.01 < 
d/ u i'c + < 100 are designated Cp and used in all calculations which in- 
clude Cp, c, Vp, and Tq. 

More than one gas constituent may have vibrational modes within the 
region 0.01 < D/uqC+ < 100. If the relaxation times differ by less than 

a factor of about three or four, the specific heats may be grouped to- 
gether with a common weighted relaxation time. This can be done because 
both 1 + s a and 1 + e-^ are slowly varying functions of D/upir + . 

Probe Shapes and Velocity Profiles 

Once the extreme values of e a and Sp are calculated, it becomes 

necessary to insert in the equations functions of the velocity profiles 
(eqs. (l6), (17), and (25)). 

A number of different probe front shapes can be used, depending 
upon the application of the pneumatic -probe pyrometer. Some of the 
possible shapes are shown In figure 6 with measured velocity profiles 
for these shapes at Mq * 0.7 (ug/uq * 1.37). The probe front shapes 

vary from flat (probes A and E) to pointed (probe C) while the constric- 
tion shape varies from a long nozzle (probe E) to an orifice (probes C 
and D). The resultant values of the relaxation function (uq/ug^Za, ob- 
tained from these profiles are shown in figure 7 as a function of the 
relaxation parameter D/uqT+. In the region D/uqT+ < 2 the relaxation 

curve depends primarily upon the exterior probe shape, while for 
D/u^t + > 4 the curve depends upon the shape of the constriction. 


nn^T. 
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The probe shape used herein is represented by probe A in figure 6. 
The subsonic velocity profiles (shown in fig. 2) were measured. Above 
sonic free- stream velocity the velocity profiles were obtained by us ing 
reference 35 and the measured subsonic profiles. The resultant values 
of (u 1 /u 2 )22 a and (u-j/ug) 2 ^ are shown in figures 8 and 9 as functions 

of the parameters u 2 /u d and D/u^t*. The values of 2-^ and T ^ are 
shown in figures 10 and 11 as functions of d/u-jT* and M^_. 


to 

i 

CQ 

O 


Use of Correction Factors to Calculate Errors 


Eather than make all the previously outlined corrections, the user 
of a pneumatic-probe pyrometer may prefer to- make certain assumptions 
as to the state of thermal equilibrium. Thus, it would be of interest 
to find the error in the calculated temperature if various assumptions 
are used, as follows: 


Method I . - Calculate relaxation times, and after finding the var- 
ious values of D/u-jT + assume all modes with D/u^t + < 100 are unre- 
laxed. Then use the formula: 


(l + e c ) 



(45) 


The error in the calculated T^ is 



» s _ + e, + 
a d c 


(46) 


Method II . - Calculate relaxation times, and after finding the 
various values of d/uj_t + assume all modes with D/vljT + > 0.01 are 

fully relaxed. Then use the formula: 


T H 




Then, the error in Tjj is given by: 



+ 


+ s d + e ( 


(47) 


“ 6 a +s b 


(48) 
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where 


h(Yp i) 

1 + Sd = ■ .- 7 ~T~ 




(49) 


Method III . - Assume all the modes are unrelaxed to the extent that 
T z i s 1 0 and that Tp f i “ Y4., then use the formula: 


f AoP 


2 r 2 


,*\2 


L III 






(50) 


The error is given by: 

^0 - T III 


T r 


*a + 6 b + 6 c + 6 d + 6 < 


where 


1 +6 <* = -hi rir 


(51) 


(52) 


Method IV. - Another method that can he used is to assume some sort 

(53) 


of "effective y" for the process. Then 

„* “(Trff ) 


Tlv ’ l T4 


where 


h(y P i) 

’ 1 + e * + 


6 b + s c + G ( 


5 

T- 


IV 


(54) 


Comparison of methods . - The general s hapes of the error curves as 
a function of a relaxation parameter D/ T \/Y'P,l®®l are shown in fig- 
ures 12 and 13. In order to calculate these curves the following assurap 
tlons were made: cp j_/R = 4.0, c-j_/R — 4.5, and c-^/R = 3.5. The param- 

eter, di ame ter divided by relaxation time multiplied by free-stream 
sonic velocity, D/ T + y / r^R^[ is used to place the curves on a common 

ordinate. A displacement of these curves upward or downward may be ex- 
pected from changes in heat capacities. 
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By using the curves in figure 12 as a guide , the following tabula- 
tion can be made of the "best” method to be used in each region of re- 
laxation. For Mj_ < 0.5: 


D/t ^/Tp <10, 

DAVrp,!^ > 10 ; 
DAVr p>1 He" < io 
D A + AAp,i Be i >10 


method I 
method II 
and (cp -j_ 
and (c-^ - 


- c 4 )/cp^]_ < 0.05, method III 
c 4 )/ci < 0.05, method III 


Above M-^ = 0.5, either D/r+^/yp or B/ U i r+ W be used to de- 

termine the best method by using the previous tabulation. 


The variation of effective y in figure 13 with Mach number should 
be noted. The way in which this term varies with Mach number depends to 
a large extent upon the degree of relaxation. 


Other Applications of Formulas 

Since 1 + e a is the mass-flow rate correction, it can be used 
to determine mg . 

To find m, as a function of Pz.l a* 1 *! Tz.Ij use equation (B25) 

“ f V ' ' 

in appendix B and equation ^17J . 

To find mg as a function of Pq and Tq, use equations (B25) , 
(33), and (37). 

Since 1 + is the pressure measurement correction, it can be 
used to determine Pg. 

To find P 2 in terms of Pz.l and Tz,l> use equations (B42) and 

(25). 

To find P 2 in terms of Po and To* use equations (B40), (25), 
(33), and (37). 
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Limitations of Equations 

The formulas in the text and in appendix B can he used for a number 
of flow processes besides those of the pneumatic-probe pyrometer. There 
are, however, certain limitations to the unmodified application of these 
formulas . 

(1) Above 4000° K the dissociation relaxation times of nitrogen and 
oxygen become short enough to be comparable to the vibrational relaxa- 
tion times. 

(2) When the mean free path of the molecules becomes comparable to 
the nozzle diameter, the continuous flow equations are no longer appli- 
cable . 

(3) Above Mj * 2.0 the rise in entropy across the shock is of such 
magnitude that the use of the expansion e^* * 1 - AS becomes doubtful. 

(4) When the equations are used to calculate long nozzle flow 
(supersonic diffusers) and M^> 2, the relaxation times cannot be con- 
sidered constant, and the evaluation must be made stepwise. 

(5) Above M-j_ * 2 the specific heats of the process cannot be con- 
sidered as simply as in equation (8) . 

For the above reasons the upper limit of usefulness of these equa- 
tions when used with pneumatic probes is given by p^ > 3X10“^ atmos- 
phere, ? 2, and 0^ ^ 3000° K. These formulas can probably be used 
with due caution for estimation purposes to pj_ = 3XiO -4= atmosphere, 

M x = 4, and Q ± = 4000° K. 


EXAMPLES 

By using the previous equations and the equations in the appendixes, 
three examples will be shown. In calculating values of 2 a and 2^ the 

measured profiles will be usedj thus, the results shown are for a probe 
with a shape similar to that shown in figure 2 (or probe A of fig. 6) 
with a 4-to-l probe-to- throat diameter ratio. 


Example 1, Products of Combustion of a Heavy Hydrocarbon and Air 

The computation of one point in this example is shown in appendix 
C. Figure 14 shows the error by using the four methods. As the free- 
stream pressure was assumed atmospheric and the temperature was assumed 
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to "be 1500° K, all dissociation effects could "be neglected. The curves 
are plotted as a function of the free-stream Mach number M-j_. The values 
of specific heats and i|r used are those given in tables I and TT . 

This example is of special interest because of the experimental data 
available at present on the use of pneumatic probes in measuring tempera- 
tures in hydrocarbon-air flames (refs. 21 to 23). The curves in figure 
14 substantiate the general conclusions for this shape of probe drawn 
from the previous figures. 

In reference 21 method III was used and was considered accurate. 

However , the probe shapes of reference 21 and this report were different. 

The ratio of probe diameters in reference 21 was as high as 8 to 1, and 
the constriction used was an orifice. Thus, compared with the present 
probe more stagnation occurred in front of the reference probe, which 
made the unrelaxed portion of the curves for longer (see fig. 7) . 

The result was that the nitrogen modes were unrelaxed, and the carbon 
dioxide and water modes were partially relaxed. For this reason the curve 
for the first three methods (applied to the probe of ref. 21) would be 
raised higher, probably by about 0.01 to 0.015, and the dependency of 
these curves upon the free-stream Mach number would be less. The curve 
for the effective y would become more horizontal and would have a value 
larger than y-p From figure 7 in reference 21 it can be assumed that the 

measured temperatures were accurate to within 2 percent. The raising of 
the error curves would place the curve for method III within this limit. 

In reference 23 an effective x which is varied from y^ at 
= 0 to y- ^ at M = 1.0 is used. The comparison between the calcu- 
lations as shown in figure 14 and the results of reference 23 should be 
direct as the sizes and shapes of the probes are similar. While the 
trend in the effective y given in figure 14 is the same as that given 
in reference 23, the magnitude of the change is different. Such a mag- 
nitude of a change in effective y in figure 14 could only be accomplished 
if the relaxation times were 100 times longer for water and carbon dioxide 
than that given in figure 5.. Thus, this discrepancy can only be attri- 
buted to other effects not considered in this paper, probably heat- 
transfer effects at low Mach numbers . 

In spite of these differences, figure 14 shows that, because of the 
high percentage of nitrogen in the combustion system, the use of any of 
the four methods does not entail any large errors over a range of subsonic 
and low supersonic Mach numbers. 


Example 2, Supersonic Flow of Combustion Products of Ammonia and Oxygen 

Figure 15 shows the errors when the various methods are used. Again, 
the errors are plotted as a function of the free-stream Mach number Mp. 
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For this case atmospheric pressure is again assumed along with a free- 
stream static temperature of 2000° K. Effect of dissociation is once 
more ignored (see fig. 6). The prohe diameter assumed was D“3/4 inch. 

In this example the vibrational modes of water could be considered relaxed 
as D/u-jT'*' >1500. The values of D/u-jT + for the nitrogen ranged from 9 
to 25. Because of the large changes in temperature and, thus, specific 
heats at these Mach numbers, the values of . fp and tyg were calculated 

from reference 32 by using a larger temperature range than that used to 
obtain table II. 

This example is of interest because it shows vhat may happen as more 
complex propellants and higher Mach numbers - are used. As the higher Mach 
numbers are reached, both methods II and I are in error. 

The downward displacement of the error in method II must also be 
noted. This displacement is due to the use of the free-stream relaxed y 
in method II where the y of the process is more nearly the y given 
by the specific heats behind the shock. 


. Example 3, Effects of Relaxation on Critical Velocity 

■ Figure 16 shows the effects of relaxation on the critical velocity 
In the first nozzle. Because of the assumption of constant critical ve- 
locity, the approximate methods cannot be used. However, the critical 
velocities were calculated for some conditions by using the iteration 
method. Figure 16 shows that the transition from the unrelaxed critical 
velocity Ujj g ^ ille critical velocity with complete thermal equilibrium 

is not a smooth one, especially in the region M]_*l and d/u^t + < 10. 


CONCLUDING REMARKS 

Equations are derived which evaluate the effects of thermal vibra- 
tional relaxation and specific-heat changes on pneumatic -probe pyrometer 
measurements In the region to 3000° K and Mach number <2.0. By using 
these derived equations, a simpler working set of equations is given 
which are shown to be adequately accurate. 

Methods are given by which the results of these equations may be 
used to determine the error In temperature measurement if different as- 
sumptions are made as to the state of thermal equilibrium of the gas in 
the probe. By using these methods, examples are given for a typical 
probe. Different assumptions must be made, depending upon the velocity, 
probe size, and gas constituents. The results are in partial agreement 
with the experimental evidence obtained in other reports. 
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The relaxation error entailed lay using any of the methods will not 
he large for the case of hydrocarbon combustion in air in the region of 
free -stream static temperature < 2500° K and free-stream Mach number < 1.2. 

The calculation of one point in one example is given as a sample 
calculation. The limitations of the equations presented and corrections 
for dissociation are discussed. The effect of partial relaxation upon 
the critical velocity is also shown. 


Lewis Plight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May 13, 1957 



APPENDIX A 


SYMBOLS 

effective nozzle area 
coefficient 

specific heat at constant pressure, (5 e/S©) 

prole outside diameter 

nozzle -throat diameter 

thermal energy (enthalpy) 

coefficient 

r±i 

n ( r ) = r ( r | 1; ) r ' 1 (e 4 . (13)) 

Mach number 

molecular weight 

mass -flow rate (eq. (6)) 

polytropic gas exponent (eq. (10)) 

total pressure 

static pressure 

specific gas constant 

entropy 

total temperature 

true total tengperature (eq. (3)) 

time 

time at point of critical flow (eq. (15)) 
velocity 
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t 

m 

o 


a 

3 

r 

r 

e 

e 

p 

2 

T 

T+ 


exponential coefficient (eq. (40)) 

correction exponent (eqs. (18), (22), (27), and (3l)) 

relaxation integral (eq. (16) ) 

ratio of specific heats, c/(c - R) 

correction factor (eq. (12)) 

static temperature 

density 

relaxation integral (eq. (15)) 
relaxation time (eqs. (l) and (2)} 


modified relaxation time, 


T 

C P,1 + C B,1 


specific-heat function, ^ ^ (eq. (9)) 

2\d0/i\ Cl / 


Subscripts: 

eff effective 

i,j mode of molecular species (used with c, t, and i|r to denote a 

particular mode or species in a general equation) 

0. 1.2, position in probe as shown in fig. 1 (used with A, c, E, M, P, 

3,4 p, T, U, if, 0, and p to denote at which physical position 

the quantity is taken) 

1, II, calculation method (used with T to denote method by which 
III, IV total temperature is calculated) 

Ihe following subscripts are used with e, 2, and T to denote type of 
correction: 

a critical flow correction 

b total -pressure measurement correction 

c correction to true total temperature 

d x correction 

e dissociation correction 
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The following subscripts are used with c } E, M, P, p, T, 9 } and t to 

denote assumed relaxation state of the variables: 

B vibrational mode when partially relaxed 

P relaxed modes when partially relaxed vibrational mode is partially 

relaxed 

R relaxed modes when partially relaxed vibrational mode is relaxed 

U relaxed modes when partially relaxed vibrational mode is unrelaxed 

X relaxed modes 

Z relaxed modes when partially relaxed vibrational mode is unrelaxed 

and specific heats are independent of temperature 

Superscripts: 

B,P, temperature state (used with E and 9 to denote temperature 

R,U state at which a mode is to be taken when that temperature state 

is not that normally assumed by the mode) 

* measured quantity (used with P, T, ft, 9 , and p to denote a 

measured quantity and/or a quantity when flow through nozzle is 
stopped) 
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APPENDIX B 

PNEUMATIC -PROBE PYROMETER EQUATIONS WITH RELAXATION 

o 

In finding the relations Between the measured quantities Pg, T£, 
and P 4 and the total temperature T 0* five equations must he satisfied: 

(1) Energy conservation (eq. (Bl)) 

(2) Critical flow (eq. (B6) ) 

(3) Entropy (eq. (B8)) 

(4) Relaxation (eq. (2)) 

(5) Conservation of mass (eq.. (B43)) 

In order to obtain a workable solution to these five equations, four 
assumptions are made : 

(1) Specific heat is a linear function of temperature: 

c i = c i,i + 2c i,i*i C e i/ e i - ^ 

where ijr < < 1. 

(2) Static pressures and temperatures are related by 

.n/n-1 


P2 f 3 4 * * * * 9 P,2\ n 


(3) Relaxation time t is constant over the process. 

(4) All modes are in equilibrium at station 1: 


0 1 0 P,1 0 B,1 0 X,1 


Energy, Velocity, and Mass Flow 

The equation for the energy at station 2 (eq. (5)) is written as 
2 u 2 

%>,2 + ®B, 2 + E X,1. + “f = E P,1 + ®B,1 + ^,1 + “^ =C P,1 T Z,1 + %,! + %,1 

(Bl) 
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or, by solving for u. 


4 ■ 2o p 


Prom assumption (2) 


0 ( ^1 *5,1 ~ ^,2 _ E P,2 ^ 

' 1 1 \ S 1 C P,1 9 1 c P,l 0 iy 


(B2) 


?2 


l -<&) 'fe) ! 


n 

. n-1 


(B3) 


where for supersonic gas flow P^ ^ is taken after the bow shock. 
By using equation (B3) the mass -flow rate is given by 


. p u = -Ssi- / flNV' 1 /£p a 2\ U_ ' L _^2_ 

l/ET Z;1 \ T Z,l7 \ e i I 


_1_ 

n-1 


m 2 = 


(B4) 


Substituting equation (B2) into (B4) and collecting terms give 
2e 


^Z.l 

m 2 = *— 




n+1 


A2?Z,1 2cp,l ( Q I \ R f ?Z,l %,!-%, 2 \ fQp.aN 11 - 1 _ /gp.zV 1 - 1 Ep,2 

V^z,i 1l R ' T z,i/ II V 0 i c p,i s i /\ 0 i/ V e i/ c p,i 0 p,; 


(2), 


(B5) 

Critical Plow and Entropy 

For critical flow dmg/dpg = 0 (eq. (6)) or because of assumption 


^/^P, 2 = 0 


(B6) 


Thus, by differentiating equation (B5) , the following equation is 
obtained: 

2 A Z,1 %,!“%, 2\ _ n+1 %*,2 _ / 9 P,2\ 2 5 ( E p,2/ c P,l 9 P,2 ^ 

e l C P,1 9 1 / n_1 C P,1 9 1 \ 9 1/ 1 S9 P,2 

29 P,2 /. 0 P,2\/ t Z, 1 • %,!"%, 2 %,2 \ dn 9 P,2 1 Se B,: 

(n-l) 2 y ®i 9 1 c P,l e l c P,l 9 lj 2 9 1 c Pi>1 50p^ ; 


(B7) 


* 
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Equation. (B7) lias three variables, Ep^gj ^B,2> n * In order to elim- 
inate n an equation for the entropy change is used. From equation (7a) 
of reference 10 or from the entropy relations of reference 6: 





+ 


dEp 



- AS 


(B8) 


Using equation (Bl) and assumption (2) gives 


1?2 n 

In — £ = m -| 
P 1 n-1 6 




By using assumption (l) , equation (B9) may he written: 


ln (s?r) E ln -f^ - <=P,1 



P 

Cx> , In -£i£ 
01 


B,2 _ 1 

1 + 2*J _ - II 


in B ? 2 

In 


-X s 

(BIO) 


Relaxation Equations 

If equations (B5) and (B7) are used to calculate the mass-flow rate 
m 2 } variables 0p g nnd 8p g must also he determined as functions 
of the state of thermal equilibrium. 

The third equation which must he satisfied is the relaxation 
equation 


dt ■% r (% “ ®b) 


( 2 ) 


If it is assumed that Og/ Cp “ Cp q/cp q over the process, equation 
(2) may he written: 
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where 


By defining 


H - 4) “ £ (% - 

(Bll) 

- t Cp ; x e - 

(B12) 

C P,1 + %1 


Er s % + ^,1 + ^,1 = E P + ®B + ^ 

(B13) 

+ %,i 

(B14) 

Er s E P + Eg + 

(B15) 

to equation (Bll) may "be written as 


<2 - %,2 + f Q [I: (E fi ' V]^ * V “ 

(B16) 


Using assumption (l) and rewriting the equation in terms of the free- 
stream Mach number and the velocity ratio u 2 /u 1 give 




^(iw M z)Bs a + ^(r^K a 

(B17) 


where 


B « 1 + 2 ( Tl ' 2 — M l) C^B “ +p) " ^p(~~ 2 — M l) ^ 
F “ (t B - t P ) - 6 *p(^- 2-" M l) C+B - 


• f* fp, , >2l 

I \A. (JL\ 

~ Jo L dt W/ -1 


U\ 2 1 (t-t f >/T + 


dt 


t P ) (B18) 

(B19) 

(B20) 
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2 


aa 



(B21) 


By using equations (B2 ) , (B7 ) , (BIO ) , and (B17) a unique solution 
for the variables 0p ^ 2 > ^ ^ or critical flow can be obtained. 

However, such a solution is a tedious process requiring double iteration 
and usually necessitating the use of an electronic computer. 


Introduction of 1 + s & 

The first step toward any simplification or approximation is to write 
the equation for mass-flow rate in such a way that the effects of approx- 
imations can bH easily seen. 

Rather than attempt to express equation (B7) in terms of (BIO), 
another form of equation (B5) can be used. The following variables are 
introduced: 


4,2 


= 2 


+ l RT z,i 


'p,i 


e 


Zj2 


L Z,1 2 
9 1 T P,1 + 1 


(B22) 

(B23) 


These would be the solutions of equations (B2) and (B7) if there 
were no partially relaxed modes (Eg = 0) and no specific-heat changes 

with temperature (Ep g - c p,i®Z, 2^* 

Inserting equations (BIO) and (2) into (B4) gives 


where 




AP. 


m 2 = 


Z,1 




i) (l + ® a ) 


-1 


*Z,2 /_X 2 


%1> 


VRT 


= Ti 


Z,l^ 


’- 1 ( r p,i + x ) 


W 1 

X -1 

r p,l 


(B24) 


(B25) 
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- 2 ASjR 


(B26) 


where 


^ - f ** a - t) 

and where , from equation (B9), 

- V &*- 

ln-gj- 

The equation equivalent to (B7) is written as 

(l + e 
d6 P,Z 


Ml + s ) 

-£L_ = o 


(B27) 


(B28) 


(B29) 


For the solution of equation (B27) it is still necessary to satisfy 
equations (B2 ) , (B7), (BIO ) , and (B17) . To simplify the solution addi- 
tional assumptions will he made. 


If the integral in equation (BIO ) , which is equal to -AS a , is 
assumed equal to zero, n and dn/dQp^g W he expressed explicitly in 
terms of and % ; 2* ®y the more simple equations which re- 

sult from the assumption of =» 0, it was found that such an assumption 
as to the entropy could he used to determine the critical velocity and 

y\q 

thus (l + e a )e a/ at Mach numbers helow = 1.4. Such an assumption 
as to the entropy leads to the rewriting of equation (B7) in the form: 



+ 


~ ^ B , 2 ®P ,2 \ 1 i" 1 / %,2 5 E B ,2 ^%,2 

c P,l e l " C p,l e l/J L 1 " R \%2 d0 P,2 + de p,2 



)]- 

(B30) 
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where from assumption (l) 


= °P,1 + %> G P,1 “ X ) 


(B31) 


SE B,2 

aa P,2 


c P.l _ c P,1 5e P,2 

°%,2 °%>,2 


(B32) 




By using equations (B30) to (B32), (B2), (B17), and (B26), the 

P/\C\ • 

value of (l + 6 )e a ‘ is found for any given state of relaxation. 

Si 

This was done using an iteration process on an IBM 650 digital computer. 
Some of the results are shown in figure 17. 

. . 34S_/R 

Also shown in the figure is the function (.1 + e a )e when this 

is computed hy the method given in the text: 


(1 + £ a ) e 


2AS a /R 


1 + 6. 




(B33) 


As is seen, this is a good approximation. 

-2AS a /R 

Next it is necessary to compute the factor e . This was 

done in the following manner: 

In the special case where ^jr = 0, AS a can he rewritten with the 
aid of the relaxation equation (2): 

' r “ (t ■ t) ■ 2 f ^ (B34) 


where 


- — f f (■ 


tf /% - %\ 2 


2 

t + 4 Lw 


dt (B35) 
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Since e " a4S a/ R « 1 . 

^ - (H 2 «if £■ 

P 7 1 


2 C 1 ° B f l r 
R a 


(B36) 


which is modified to the following when there is a change in. specific 
heat with temperature: 


-2AS a /R 


* 1 


- (H 2 M 1 f ^ + 2 (H 2 M ?) C*B - tp>] (BS7) 


In most cases equation (B36) is sufficiently accurate. 
Combining equations (B33) and (B36) gives 


6 a = ®a 


’ u + ( x ' (vv 2 - 1) (Sa ’ E ' 6a ’ D> ‘ ^ r ' 


and repeating equation (B24) results in 


a 

(B38) 


APp 1 - - i ■ - . . 

■“2 ° VhCr>,l)(l + «a) 


(B24) 


Values of (u 1 /u 2 ) 2 S a as a function of d/u^t* and u^/u^ are 
shown in figures 7 and 8. Values of (u-^/ug)^ r a as a function of 
d/uj_t + and ug/u^ are shown in figure 9. 


Introduction of 1 + e-^ 


The flow in the prohe was stopped, and P 2 was measured. This 
measured value of P 2 is used in the equation for mass-flow rate (eq. 
(B25)) which is now written: 


m2 


lfe>,l)(l + £ a)Cl + ®h) 


(B39) 
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o 

OJ 


where 


1 + £ b = 



(B40) 


or, by using the entropy equation (B9) and the form of (B26) and (B28) 


where 


1 + *, 


■ft 1 ) 


s 2cp A /E , e .2fi* 2£Sh/R 


xp> 


2/ 




C P l^P f^~2 


-) 


+ ^i 


0 * 

m JS*2 

0, 


P.1 L 

m-g- L 


1 + 1(f. 


'n* 

% 

B \ 6-. 


(B41) 


(B42) 


(B43) 


Since it is not necessary to determine ug/u-^ for this case, *2 
may be obtained directly from a modified form of equation (B17) in which: 

(l) %^2 ^- s r epla ce< 3- "by T* 


(2) 2 a is replaced by 2-^ 

(3) 2^ is replaced by 2^ 

By using the equation for conservation of energy 3- e 'k er ~ 

mined. Thus, T|, fl* 2 , and P* can be calculated. 

The correction eq.ua ti on (25) may be obtained in the following 


manner: 
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For the special case of if = 0 the pressure Pg is obtained from 
the entropy equation (B9) : 


In 


Pg f /l 1 \ f d0g 

si “A ° b (% ■ %) + A ° e + 0p 

“X ° B i0B ft " ft) 


d0p 


9 B,2 6 P,2 

+ C B 121 + C P ln 


(B44) 


(B45) 


Using once more the relaxation equation (25) as was done in equation (B34) 
and expanding the last two terms in the previous equation give 



, *1 , V 

In — * c In -r- + Oq 
Pq a 





(B4,7) 


where 



(B48) 


^ 2 
The pressure Pg/p-^ is given In the form A + B (Z-^ + F^) which, 

when inserted into equation (B40) and in turn normalized for the ex- 
tremes, is the same as equation (25): 


6 b = 6 b,U + (1 " z b " r b^ ^ e b,R " 6 b,u) 


(25) 


f Thp calculated values of are shown in figure 11 as a function 

of d/u]_t + and M-j_. In the range 0.3 < M^_ < 1.4 the values of 
for d if ferent specific heats were equal to within 10 percent of those 
shown. 
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Figure 10 shows the calculated values of 2-^ as a function of 
D/u-jTf*". A comparison of the function (l + calculated from 

equation (B4l) and from equation (25) is shown in figure 18. Ihe agree- 
ment is good. 


Flow through Second Nozzle 

The second nozzle is used as a method of measuring the mass-flow- 
rate. If it is assumed that the gas temperature at this nozzle is lcw 
(e 3 < 500° K) , the usual equations for critical flow through a nozzle 

can he used (refs. 19 to 23): 


m 


a 4 f* 


4 


'RTJ 


k(rj 


(B49) 


Final Equation 

If mass is conserved through the probe, 

m 2 = m 4 

Equations (B39) and (B42) are combined to yield 


or 


Tz ' 1 '(®^Si Cl + ea)(1 + s ’ j) 

T ° = (SfJ ^ (1 + 6 * Kl + e > )(1 + 


where 


1 + s„ = 


c 


Z,1 


( 8 ) 


(B50) 


( 12 ) 


(B51) 


Above * 1.4 the entropy change AS a can no longer be treated 
as a correction to the function 1 + e a because this entropy rise will 
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change the ratio ug/u^ as found from equations assuming AS a = 0. The 

velocity ratios may .he found hy using as end conditions the velocity 
ratios calculated hy the formulas of reference 32 for constant specific 
heats: 


u r 


u 


T 2 ”! + 1 


r + i 


(B52) 


M? 


This was done for the points calculated in figure 17, for = 2.0, 

and for the supersonic velocity ratios shown in figure 16. As is seen in 
figure 17 the approximation of 1 + * a /e - ^ a given by eq.ua tion (B32) is 

still good. 

To find another simplifying assumption can he made for 

Mp > 1.4. In this region almost all the entropy rise is due to the bow 
shock, and, thus, r a can be calculated for the normal shock alone. -Also, 
above = 1.4 and when D/upT + > 1.0, AS a “ Z£>p. For these reasons, 
and e-^ may be combined to give, for Mp > 1.4: 

6 a + 6 b 6 a,U + f 1 ' (6a > R “ 6& > u) + 6b ' U + 


(*" - ^b R b,subsonic^ 6 h,R “ 6 b,U^ 


(B53) 


This approximate formula seems- to be good in the region 1.4 < Mp < 2.0. 
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APPENDIX C 

SAMPLE CALCULATION FOR FIGURE 14 

Given 0-^ = 1500° K, P^ = 1 atmosphere, M^. ^ = 0.7, gas constituents 

of 13.6 percent carbon dioxide, 13.8 percent wa£er, and 72.6 percent 
nitrogen, a probe diameter of 1.27 centimeters (l/2 in.), and 0g < 500° K. 

First, it is necessary to estimate the value of D/u-jT+ of al 1 the 
gas constituents so as to decide which vibrational modes are to be con- 
sidered relaxed, partially relaxed, and unrelaxed. 

Adding the specific heat by using the percentage composition and the 
values given in table I gives 

c^R = 4.760, r L = 1.2650, c^/R => 3.797, and y~ 4 “1.386 

= 7.4X10^ cm/sec 

Dividing the relaxation times in figure 5 by the proper partial pressures 
gives: 

t N -N = 340x10”® sec (using eq. (40)) 

2 2 

t N -HO* 36x10"® sec (using eq.. (4l) ) 

2 2 

“ 33x10"® sec (using eq. (42)) 

a 2 

X C0 -CO * °* 29xl0_6 sec (using eq. (40)) 

2 2 

T C0 2 -H 2 0 “ °* 14x10-6 see (using eq. (41)) 

.*. t s 0.09x10"® sec (using eq. (42)) 

CO 2 

t H 2 0 “ 0.05x10”® sec (us ing eq. (40)) 

(D/u-[T + ) N2 » 0.5 

(d/u-jtt + )co 2 “ 19 ° 

(d/u 1 t+) H2 q “340 
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In the region 0.1 <d/u2_t + < 100 a mode is partially relaxed. ' Thus, 

only nitrogen will he considered partially relaxed, while water and carbon 
dioxide will be considered fully relaxed. 

Using tables I and II or equation (8) and reference 24 gives 

C-Tj ■* C-r> *1 Ct 


T,1 . ___ °B,1 U X, 1 

-g- = 4.258, — ^-is0.502, — 


0 


^- 0 . 078 , - 0 . 122 , _ 0 


T. 


Z 1 

■ =* 1.0760 (using eq. (27)) 


%2 


“0.9342 (using eq. (19)) 


( u Z,2/ u l) 2 * 1 * 94:5 (using eq. (14)) 
Pjj a -0.00513 (using eq. (20)) 
0 U,2/ 0 l“ 0,954tl (using eq. (2l)) 

1 + “1.0001 (using eq. (19)) 

Pp = 0.490 (using eq. (23)) 

0p 2/ 8 l“ 0,94=09 (using eq. (24)) 

1 = e a ^p = 0.9785' (using eq. (22)) 
P* “0.0059 (using eq. (28)) 

%,2/ 9 l “ 1*0759 (using eq. (29)) 

1 + s-jj y “1.0000 (using eq. (27)) 

Pp “0.516 (using eq. (3l)) 
T R, 2 / S l“ 1,0678 (using eq. (32)) 

1 + e-jj pa 1.0015 (using eq. (30)) 
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Using the estimated value of d/u^t* gives 

M^V 0,62 C usin s fi e* ? ) 

(u^/ug)^ r a “0.06 (using fig. 8) 
2^ “0.04 (using fig. 9) 

“0.27 (using fig. 10) 

6 “ 0.006 (using eq. (17)) 

a 

e^s: 0.001 (using eq.. (25)) 


For method I: 


For method II: 



“0.007 


(using eq.. (45)) 



“0.023 


(using eqs. (33), (48) , and (49)) 


For method III: 


T - T 

— -- — “ -0.025 (using eqs. (33), (5l), and (52)) 

■*■0 

For method IV a graph of y plotted against h(y) or the graph in figure 
11 of reference 23 (where h(y) s Fg(y)2) is used. 

y eflf “1.307 (using eqs. (13), (33), and (54)) 
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TABLE I. - TRANSLATIONAL PLUS ROTATIONAL AND VIBRATIONAL 
VALUES OF SPECIFIC HEATS FOR DIFFERENT PURE 
GASES (IDEAL GAS VALUES) 


[Data obtained from ref. 24.} 


Mode . 

Tempera- 

ture, 

°K 

c/R 

(a) 

Carbon 

dioxide 

Water 

Nitro- 

gen 

Oxygen 

Carbon 

mon- 

oxide 

Translational 

plus 

rotational 

All 

tempera- 

tures 

3.5000 

4.0000 

3.5000 

3.5000 

3.5000 

Vibrational 

500 

1.8671 

0.2345 

0.0578 

0.2396 

0.083 

Vibrational 

1000 

3.0318 

0.9569 

0.4326 

0.6958 

0.491 

Vibrational 

1500 

3.521 

1.6526 

0.6919 

0.8975 

0.736 

Vibrational 

2000 

3.758 

2.1460 

0.8268 

1.0436 

0.859 

Vibrational 

2500 

3.893 

2.4727 

0.9047 

1.1808 

0.929 

Vibrational 

3000 

3.984 

2.6945 

0.9545 

1.3062 

0.976 


Temperature , c/R for air 

OK (a) 



All trans- 
lational 
and 

rotational 

modes 

All vi- 
brational 
modes 

All modes 
except nitro- 
gen vibration- 
al modes 

Nitrogen 

vibra- 

tional 

modes 

500 

3.4906 

0.0959 

3.5414 

0.0451 

1000 

3.4906 

0.5044 

3.6374 

0.3376 

1500 

3.4906 

0.7287 

3.6802 

0.5391 

2000 

3.4906 

0.8855 

3.7109 

0. 6452 

2500 

3.4906 

0.9550 

3.7397 

0.7059 

3000 

3.4906 

1.0203 

3.7660 

0.7448 
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TABLE II. - VALUES OP SPECIFIC-HEAT FUNCTION (IDEAL GAS) 


Mode 

Temp era - 
ature, 

°K 

(al 

Carbon 

dioxide 

Water 

Nitro- 

gen 

Oxygen 

Carbon 

mon- 

oxide 

Translational 

plus 

rotational 

All 

tempera- 

tures 

0 

0 

0 

0 

0 

Vibrational 

1500 

0.140 

0.544 

0.402 

0.328 

0.346 

Vibrational 

2000 

0.093 

0.299 

0.243 

0.270 

0.209 

Vibrational 

2500 

0.069 

0.270 

0.169 

0.281 

0.151 

Vibrational 

3000 

0.058 

0.206 

0.129 

0.267 

0.123 

Temperature , 
°K 

V for air 

(a) 


All trans- 
lational 
and 

rotational 

modes 


1500 

2000 

2500 

3000 


All vi- 
brational 
modes 


0.368 

0.298 

0.198 

0.193 


All modes 
except nitro- 
gen vibration- 
al inodes 


0.020 


0.020 

0.022 

0.024 


Nitrogen 

vibra- 

tional 

modes 


0.402 

0.243 

0.169 


0.129 


a 


1 

2 






Adiabatic 
6 1 agnation. 


0 


Tree Critical 

stream flow 



First 

nozzle 


Figure 1. - Schematic diagram of 
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Velocity ratio, u/u-^ 
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Figure 2. - Probe shape and experimentally obtained velocity profiles. 
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0 (Station l) Time, t — »- 



(Station 2) 


Velocity, u 



Vibrational 

tengjerature 


Relaxed, 0 R (i«t f ) 


— Partially relaxed, 

Unrelaxed, 6-^ (t» 



Temperature of 
relaxed modes with 
vibrational mode — 

Relaxed, Sg 

Partially relaxed, Up 

Unrelsxed, 


Figure 3. - Schematic diagram showing effects of relaxation on temperature. 




Change in gas constant, (R^ - R^/E^ = 
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Temperature, °K 

Figure 4. - Effect of dissociation on variation of gas con- 
stant with, temperature for three combustion processes at 
atmospheric pressure. 
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Relaxation time t for 1-atm partial pressure of molecule j, microsec 




Relaxation time t for 1-atra partial pressure of molecule J, microsec 
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.04 .06 .08 .10 J, 12 .14 .16 .18 

Temperature, (®K )V3 

I I I I I 

5000 2000 1000 500 273 

Temperature, °K 

(b) Chlorine, carbon dioxide, nitrous oxide, and water. 

Figure 5. - Concluded. Experimental relaxation times as functions of 
temperature . 
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Figure 7. - Variation of relaxation function Z a with relaxation parameter far five probe shapes 
shown in figure 6. Free- stream Mach number, * 0.7; velocity ratio, 1.57. 
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Figure 8. - Variation of relaxation funotion with velocity ratio and relaxation parameter. ^ 

Typical probe shape (probe A of fig. 6 or probe of fig. 2). W 
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.001 ... .01 .1 1 10 

Relaxation parameter, D/u]_t + 


■ Figure 11. - Variation of relaxation function with, free-stream Mach number and relaxation 
parameter. Typloal probe shape (probe A of fig. 6 or probe shape of fig. 2); 

0.3 < Mach number <1.4. 
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(c) Mach number, *1,4. 

Figure 12. - Concluded. Variation of temperature error with relaxation parameter for three methods of 
calculation, y^ ® 1.333: - 1.285j « 1.400; f p “ ” 0- 
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